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Abstract. We develop a new general algorithm for finding a regular tight-binding lattice 
Hamiltonian in infinite dimensions for an arbitrary given shape of the density of states 
(DOS). The availability of such an algorithm is essential for the investigation of broken- 
symmetry phases of interacting electron systems and for the computation of transport 
properties within the dynamical mean-field theory (DMFT). The algorithm enables us 
to calculate the optical conductivity fully consistently on a regular lattice, e.g., for the 
semi-elliptical (Bethe) DOS. We discuss the relevant /-sum rule and present numerical 
results obtained using quantum Monte Carlo techniques. 



1. Introduction 

Microscopic studies of strongly correlated electron systems require methods 
which take the Coulomb repulsion between electrons explicitly into ac- 
count. Many features of such systems can be modeled using the single-band 
Hubbard model [1] 



H = H + Hint = % ^Ria^a + U Y1 "Hrf"**! ' W 



where the operators c l R , a and c R , a create and destroy electrons of spin 
a on site Ri, respectively; fij^a- measures the corresponding occupancy. A 
general nonperturbative treatment of this model is only possible in the limit 
of infinite dimensionality where the dynamical mean-field theory becomes 
exact: due to a local self-energy the model reduces for d — ► oo to a single 
impurity Anderson model plus a self-consistency equation. For homoge- 
neous phases, local properties then only depend on the lattice via the 
noninteracting DOS p(e) := jj J2k ~~ e fc)- I n the simplest case of uniform 
nearest-neighbor hopping on a hypercubic (he) lattice, the dispersion reads 
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Figure 1. Bethe lattice: a) conventional tree level picture for coordination Z = 4, b) 
stacked Bethe lattice, c) redefined Bethe lattice as a regular lattice with cubic symmetry: 
here, symbols denote equal hopping matrix elements relative to a fixed site (central circle). 



e£ c = — 2t J2a=i cos(k a ) which for the proper scaling t = t* /y/2d leads to a 
Gaussian DOS p hc (e) = exp[-e 2 /(2t* 2 )]/(^t*). 

Vertex corrections to the optical conductivity <t(cj) vanish [2] in the 
limit d — > oo so that it may be expressed in the isotropic case as [3] 

a(u) = a o rdep(e) [°° du' A £ (u>)A e (u> + u f M ~ ^ + ^ . (2) 

J —oo J —oo OJ 

Here, A e (uj) is the "momentum" dependent spectral function, nf(a/) is the 
Fermi function, oo = ^p-y 3 and 

~p^) ■= jj E l^l 2 ^( e - ^) =■ (l^l 2 )( £ ) p^) ■ (3) 

Note that the frequency- and interaction-dependent part in (2) depends on 
the lattice only via the DOS while the explicitly lattice-dependent part p{e) 
is universal, i.e., independent of interaction U, filling, and temperature. In 
the hypercubic case, the momentum dependence of the Fermi velocity 
becomes irrelevant in (3) for d — > oo; for unit hopping and lattice spacing, 
one observes /0 hc (e) = p hc (e)- As a consequence, the optical /-sum is then 
proportional to the kinetic energy: 

l o °°du;a(u) = -^(Ho}. (4) 

However, the he DOS is unbounded which is hardly compatible with the 
single-band assumption. In fact, no regular lattice model with sharp band 
edges in d — » oo could be constructed so far. In this situation, many DMFT 
studies have focussed on the so-called Bethe lattice which is not a regular 
lattice, but a tree in the sense of graph theory as shown in Fig. la. The 
semi-elliptic DOS p(e) = y/A — e 2 /(2ir) of this model (for Z — > oo) fixes 
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the local properties of the model; transport, however, is a priori undefined. 
A derivation of <j(u>) directly for the Bethe tree (using the level-picture 
Fig. la) by Chung and Freericks [4] is still incomplete [5]; up to a factor of 
3, the same expression p(e) oc (4 — e 2 )p(e) was obtained [6] in a heuristic 
scheme by enforcing the he /-sum rule (4). An alternative direct approach 
[7] fails to describe the coherent transport expected in the metallic regime. 

The local DMFT problem is unchanged when a finite number of hopping 
bonds per site are added. Therefore, the periodically stacked Bethe lattice 
(Fig. lb) is still a Bethe lattice in the DMFT sense; potentially coherent 
transport is then well-defined (only) in stacking direction [8]. A semi-elliptic 
DOS also results from fully disordered hopping on lattices of arbitrary 
topology; in this case, <t(uj) is incoherent [9]. 

We will in the following construct and evaluate a new definition for 
ct(uj) compatible with a semi-elliptic DOS. This definition is unique by 
being based on a regular lattice as illustrated in Fig. lc and by leading to 
an isotropic conductivity which is coherent in the noninteracting limit. 

2. General Dispersion Method 

We rewrite the translation-invariant noninteracting Hamiltonian, 

^ = 11^4/^ = 1^)^' ( 5 ) 

where contributions to the dispersion may be classified by the taxi-cab 
hopping distance ||r|| = J2a=i \ T a\- 

oo 

e(fc) = £ e D (fc), e D (k)= t ^' h - ( 6 ) 

D=l \\t\\=D 

In high dimensions, only vectors of the form r = Y^f=i e a.i with pairwise 
different directions CKj need to be considered. This follows from the fact 
that the fraction of neglected vectors (with |r ■ e a \ > 1 for some direction 
a) vanishes as 1/d. Furthermore, the considered vectors are of minimal 
Euclidean length |r| hinting at maximal overlap, i.e., largest \t T \ for fixed 
taxi-cab distance ||r|| and fixed D. 

By deriving a recursion relation for eo(k) we have established that 

OO J.* 

D=l v L> - 

Using the orthogonality of the Hermite polynomials, one may express the 
hopping matrix elements in terms of the transformation function J~{x) : 

t*D = -y^y j°°JeHe)Ke D (e)e-^l\ (8) 
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Specializing on the case of a monotonic transformation function J-(x) (with 
derivative ^'(x)), we can write 



which leads to 



JF-!( e ) = y/2 err 1 ^2 J' de'p(e') - lj . (10) 
Furthermore, the Fermi velocity = Ve^ can be computed: 

„ - T' ( < e))v hc - P hC (^err 1 (2/- 00 d e y g Q-l)) 

v k -J-{J- {e))v k — v k . (11) 

A practical application of the general formalism proceeds as follows: 

1. compute .F -1 (e) from arbitrary target DOS p(e) using (10) 

2. invert function (numerically or analytically) to obtain J-(e) 

3. evaluate transport properties, e.g., (\v k \ 2 )(e) or /3(e) using (11) 

4. optionally determine microscopic model parameters t* D using (8) 

The only choice inherent in this procedure beyond the usual assumptions 
for large dimensions is contained in step 1 which by construction produces 
a monotonic transformation function T . The optical /-sum rule reads 

/foW^u,) = g = g i\T"{T-\,)) - T-\e)T\T-\,)) 

(12) 

Here, the first equality follows from (2), i.e., is generally valid within the 
DMFT while the second expression in terms of the transformation is specific 
to the formalism developed within this section. 

Exam-pit: Flat-band System One interesting limiting case of a monotonic 
transformation function which can be treated analytically is the step func- 
tion F(x) = 20 (x) — 1 corresponding to a flat band DOS of the form 
p(e) = {S(e— 1) +<5(e+l))/2. For this case, the hopping matrix elements 
read (trivially, t\ n = 0): 

*$n+i = (-Wf feS ™ (-i) n (™r 3/4 - (13) 

V 7r v /(2n + 1)! 

The asymptotic exponent —3/4 is only slightly smaller than the threshold 
value —1/2 required for a finite variance f^° 00 dee 2 p(e) = J2d^*d 2 - F° r a 
rectangular model DOS, t* D ~ 2 -n n~ 3 / 4 already decays exponentially fast. 
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Figure 2. a) The average squared Fermi velocity (|ufc| 2 )(e) is constant for the hypercubic 
lattice (or for the ^-component of a stacked lattice); in contrast it vanishes for the isotropic 
lattice defined in this work. For comparison, the form suggested by Millis is also shown, 
b) Resulting function p(e) of the full isotropic model (solid line) in comparison with 
truncated models (D ma x = 3 or D ma x = 5), evaluated in finite dimensions 5 < d < 100. 



3. Application to the "Bethe" semi-elliptic DOS 

In this section, we will apply the new formalism to the Bethe semi-elliptic 
DOS in order to determine a corresponding tight-binding Hamiltonian de- 
fined on the hypercubic lattice with the same local properties as the Bethe 
lattice (with NN hopping) in the limit d — > oo. From (11), we derive the 
average squared Fermi velocity defined in (3) in closed form: 



I l2\ t \ 27T 

>( e ) = i — ^ ex P 



_ 2 Lf-i ^Vl ~ g V 4 + 2 arcsin(6/2) j\ 2 



(14) 



Here, we have used the fact that (|i>fc| 2 )(e) is effectively constant (and equals 
1 for unit variance and lattice spacing) in the hypercubic case. The result 
(solid line in Fig. 2a) has all the qualitative features expected for this 
observable in any finite dimension: (|i?fc| 2 )(e) is maximal near the band 
center, strongly reduced for large (absolute) energies and vanishes at the 
band edges: states at a (noninteracting) band edge do not contribute to 
transport. The violation of this principle in the stacked case (dashed lines 
in Fig. 2a), which corresponds to an application of the he formalism to 
the Bethe DOS with (|i>fc| 2 ) constant up to the band edges, is clearly 
pathological. Therefore, our method has not only the merit of yielding 
isotropic transport, but also of avoiding unphysical behavior. 

In order to determine the microscopic model, we have to apply (8) to 
the numerically evaluated transformation function T. Again, the scaled 
hopping matrix elements fall off exponentially fast: only a fraction 10~ 3 
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Figure 3. Numerical results for the half-filled Hubbard model with semi-elliptic DOS for 
d — > oo in the paramagnetic phase at T — 0.05. a) Local spectral function A(cu) obtained 
from QMC (using a discretization At — 0.1) and MEM. b) Optical conductivity u(u>) for 
the isotropic "redefined Bethe lattice". The inset shows the partial /-sum J Q ' Auj' 'a(uj'). 



of the total energy variance arises from hopping amplitudes beyond third 
nearest neighbors and only a fraction 10 -6 results from hopping beyond 9 th - 
nearest neighbors. This result suggests that properties of the model should 
be robust with respect to truncation. In fact, /5(e) (and consequently the 
definition of (j{u))) hardly changes when hopping is cut off beyond 3 rd or 
5 th nearest neighbors, even when evaluated in finite dimensions as seen 
in Fig. 2b. This behavior is very general so that results for cr(w) of a 
local theory in finite dimensions will depend on d predominantly via the 
interacting DOS A(uj) and only very little via p(e). 

The local spectral functions for T = 0.05, i.e., slightly below the critical 
temperature T* ~ 0.055, are shown in Fig. 3a as obtained from QMC/MEM 
[5]. In the metallic phase, the spectral density at the Fermi level (to = 0) 
is approximately pinned at the noninteracting value p(0) = 1/ir ~ 0.32 
for U < 4.4. The quasiparticle weight decreases drastically and a shoulder 
develops for U > 4.6 before a gap opens for U > 4.8. An application of 
(2) to these spectra for the isotropic model characterized by (14) yields 
the estimates for the optical conductivity a(u) shown in Fig. 3b. A low- 
frequency Drude peak (of Lorentzian form) and a mid-infrared peak at 
ijj ~ U/2 are present in the metallic phase and decay towards the metal- 
insulator transition at U ~ 4.7. For large U, the optical spectral weight 
concentrates in incoherent peaks at uj ~ U. The inset of Fig. 3b shows the 
partial optical /-sums. As expected, both the contribution of the Drude 
peak and the total /-sum decrease for increasing U. 

Figure 4a shows the impact of the definition for (j{uj) on the results 
for U = 4.0. Our isotropic model yields by far the largest contributions 
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Figure 4- a ) Optical conductivity <j(uj) for T — 0.05 and U = 4.0 for the new 
isotropic model with semi-elliptic DOS in comparison with disordered and stacked models 
consistent with the same DOS. b) Optical /-sum for T — 0.05. 

at small uj. While the stacked model leads to otherwise similar results, 
the low-frequency form of ct(lo) is qualitatively different in the disordered 
case (where a Drude peak is absent even for U — > 0). As seen in Fig. 4b, 
the /-sum is generally larger for the isotropic than for the stacked Bethe 
lattice, in particular in the metallic phase (while the disordered case is 
in-between). These differences can be attributed to the enhanced squared 
Fermi velocity in the isotropic model: The enhancement is largest near the 
Fermi surface, at e = by a factor of tt/2 ~ 1.57. A corresponding increase 
is expected of the Drude peak for small enough U and T when transport 
is dominated by states with e « 0. Since energy eigenstates spread out in 
momentum space at large U, the enhancement reduces (in general) to the 
integral j^° 00 dep(e) = J2d=i Dt* D 2 which evaluates here to 1.05406. 

In all three cases, the proportionality (4) of the /-sum to the kinetic 
energy as characteristic for the he lattice is clearly violated. This is true even 
for the stacked case (which is otherwise similar to the he case): while the 
/-sum is here proportional to the contribution to the kinetic energy associ- 
ated with hopping in current direction [8], / °°da; a xx (u>) = — o"o(T x )/4, this 
contribution (which is negligible in the limit Z — > oo) is not proportional 
to the total kinetic energy in this anisotropic case. A more relevant sum 
rule is derived from (12): J^dua^u) = ^r^ijE^)- 

4. Conclusion 

We have presented a new general method for constructing regular lat- 
tice models with hypercubic (he) symmetry, i.e. isotropic optical transport 
properties, in large dimensions. Previously, calculations of the optical con- 
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ductivity a{oS) of the Hubbard model in the limit d — > oo had been restricted 
to the hypercubic lattice (using NCA [3] or QMC [10]) or have ignored the 
lattice dependence: Applying the he formalism to the Bethe semi-elliptic 
DOS, Rozenberg et. al [11] overlooked violations of the he /-sum rule. 
All previous approaches specific to the Bethe DOS were associated with 
anisotropic or incoherent transport; the most interesting of these [4, 6] 
have not yet been linked rigorously to microscopic models. 

Our method yields the first derivation for cr(u) consistent with a semi- 
elliptic DOS that implies isotropic transport which is fully coherent in the 
noninteracting limit. This reinterpretation of the "Bethe lattice" (in the 
DMFT sense) as an isotropic, regular and clean lattice and the demon- 
stration that the associated transport properties are robust (with respect 
to finite dimensionality or hopping range) removes, finally, the pathologies 
previously associated with the DMFT treatment of transport in connection 
with non-Gaussian DOSs. At essentially no additional cost, the method 
can also be used for computing properties such as transverse conductivities 
and thermopower; these vanish, however, in the particle-hole symmetric 
case considered in this paper. Our numerical results have shown that the 
precise definition of <r(u>) does matter, in particular within the metallic 
phase where transport is potentially most coherent. We have also found a 
general DMFT expression for the /-sum rule as well as a form specific to 
our new approach. 
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